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We investigate the dark matter annihilation in the brane-world and quintessence scenarios, in which the modi- 
fied cosmological expansion rate can enhance the thermal relic density of dark matter. According to the observed 
dark matter abundance, we constrain the thermally averaged annihilation cross section (crv) in these two sce- 
narios. In addition, the big bang nucleosynthesis and the partial-wave unitarity are also used to place bounds on 
(av). It is found that both scenarios can lead to a large annihilation cross section, so they can be used to explain 
the recent PAMELA, ATIC and PPB-BETS anomalies. 

PACS numbers: 95.35.+d, 98.70.Sa 



I. INTRODUCTION 

The existence of dark matter is by now well confirmed [Ql]. 
The recent cosmological observations have helped to establish 
the concordance cosmological model where the present Uni- 
verse consists of about 73% dark energy, 23% dark matter, 
and 4% atoms 101 ■ However, in the standard model of particle 
physics, there is no candidate for dark matter. Understanding 
the nature of dark matter is one of the most challenging prob- 
lems in particle physics and cosmology. One particularly at- 
tractive class of dark matter candidates is provided by Weakly 
Interacting Massive Particles (WIMPs). So far, we still do not 
know the WIMP's mass m and the thermally averaged annihi- 
lation cross section {crv). 

In the standard cosmology, the observed dark matter abun- 
dance Qo/j^ = 0.1131 ± 0.0034 |:2:] requires (cry) 3 x 
10"^^ cm^ sec"', if the i-wave annihilation is dominant lU. 
In this case, (crv) is insensitive to the temperature T of the 
Universe. In fact, the thermal relic density of dark matter de- 
pends not only on their annihilation cross section {crv), but 
also on the cosmological Hubble expansion rate during the 
era of dark matter production and annihilation. So, any de- 
viation from the standard cosmology, at this epoch, will lead 
to some variation of the thermal relic density of dark mat- 
ter. Actually, before the big bang nucleosynthesis (BBN), the 
Hubble parameter remains unclear. It is likely that some un- 
known mechanisms govern the evolution of the Universe in 
the epoch. For example, the brane-world scenario 101 and the 
quintessence scenario with a kination phase Jstl, may play a 
significant role in the pre-BBN era. In these two scenarios, 
the enhancement of the Hubble expansion rate will lead to the 
enhancement of the thermal relic density of dark matter |6il . 

For the observed relic density, the modified Hubble param- 
eter scenarios possess a larger annihilation cross section than 
(crv) R! 3 X 10"^^ cm^ sec"' in the standard cosmology. There- 
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fore, many dark matter models can derive the larger param- 
eter space when we choose different Hubble parameter. On 
the other hand, some models, which are excluded in the stan- 
dard cosmology, can work well in these scenarios. It is very 
important for us to derive bounds on the thermally averaged 
annihilation cross section {crv). In this paper, we try to give 
a comprehensive constraint on the dark matter annihilation 
cross section in the brane-world and quintessence scenarios. 
Besides, the BBN and partial-wave unitarity can also be used 
to constrain the parameter space of the above two scenarios. 
This paper is organized as follows; In Sec. II, we investigate 
the constraints on (crv) from the relic density, BBN and uni- 
tarity bounds. We also discuss the recent PAMELA, ATIC and 
PPB-BETS anomalies. Some discussions and conclusions are 
given in Sec. III. 



II. CONSTRAINTS ON DARK MATTER ANNIHILATION 
CROSS SECTION 

The evolution of dark matter abundance is given by the fol- 
lowing Boltzmann equation : 



dY s(x) T T 

';{(TV){Y^ - YIq) 



(1) 



dx xH 

where x = m/T, Y = n/s{x) denotes the dark matter number 
density, and Yeq is the equilibrium (EQ) number density. 



-r-7—x K2(x) 
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with K2(x) the modified Bessel functions, and the total 
number of effective relativistic degrees of freedom. For the 
internal degrees of freedom of dark matter particle, we take 
gi = I. The entropy density s(x) is given by 



s(x) 



(3) 



45 x3 ■ 

In the standard cosmology (SC), the Hubble expansion rate H 
is written as 
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where Mpi = 1.22 x lO'"* GeV is the Planck (PI) mass. Using 
the result of Eq. ([1]), Yq, one can obtain the dark matter relic 
density Qb/z^ 10]: 



Q.Dh^ = 2.74 X 10 



8 



GeV 



(5) 



The thermal average of annihilation cross section times the 
"relative velocity", {crv), is a key quantity in the determination 
of the cosmic relic abundances of dark matter. In the standard 
cosmology, one can use approximate formulas to calculate the 
dark matter relic density. Then, (crv) a; 3xlO"^^cm^ sec"' can 
be obtained for the ^-wave annihilation lH,!^. If the Hubble 
expansion rate H deviates from that of the standard cosmol- 
ogy, one should numerically resolve the Boltzmann equation 
in Eq. ([T]i- In this case, we have three parameters: m, {crv} and 
H. For the WIMP, its mass m should be roughly between 10 
GeV and a few TeV for annihilation cross section of approxi- 
mately weak strength. So, we take 10 GeV < m < 10 TeV. In 
this paper, we only focus on the i-wave annihilation of dark 
matter. Thus, (crv) is assumed to be a constant. In terms of 
the partial-wave unitarity, one can derive the unitarity bound 
on (crv) H]: 



(o-v) < 



3 X 10 ^^cm-' sec 
(OT/TeV)2 



(6) 



We shall show that this limit will exclude part of the parameter 
space. Using the observed dark matter abundance Q/j/i^ = 
0.1131 + 0.0034 Q], we discuss the constraints on the above 
three unknown parameters in the following subsections. 



A. The brane-world scenario 

In this subsection, we consider the dark matter annihilation 
in a brane-world scenario, in which the standard model parti- 
cles are assumed to be confined on a 3-brane. In the early Uni- 
verse, extra-dimension effects might play an important role. It 
is of interest to study how the physics of extra dimensions may 
affect the dark matter annihilation. 

Focusing on the case with one extra dimension compact- 
ified on a circle, the effective four dimensional Friedmann 
equation is 



Stt 



3M; 



-P 



1 + ^ 

Pc 



(7) 



where p - n-g^T'^ /3Q is the usual energy density of radiation, 
and Pc = 2A is the critical density, with A the brane tension. 
The critical density pc can be expressed as 



Pc = 



967rM^ 
Ml 



(8) 



where M5 denotes the true gravity scale of the five dimen- 
sional theory. Here, we have ignored the dark-radiation term 
in Eq. Q, considering the BBN constraint. On the other hand, 
the BBN can also constrain the critical density pc, 
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FIG. 1: The allowed region of (crv) and m for different M5 in the 
brane-world cosmology from the observed relic density Sioh^ = 
0.1131 ± 0.0034. Here, we also consider the unitarity bound in Eq. 
^ and the BBN bound in Eq. dlOl l. The region with dotted bound- 
ary denotes the allowed range of (crv) in Eq. i20\ from the PAMELA 
experiment. 



Then, we straightforwardly derive the BBN bound, 



Ms > 1.1 X 10'*GeV 



(10) 



It should be mentioned that the precise measurements of the 
gravitational law in submillimeter scale give a more strict limit 
M5 > 1.1 X 10^ GeV II^l. However, this constraint is model 
dependent ifllll . So we only consider the BBN bound M5 > 
1.1 X 10"* GeV in the following analysis. 

Using the modified Friedmann equation Eq. O and the 
observed relic density Q/j/i^ - 0.1131 + 0.0034, we numer- 
ically solve the Boltzmann equation. Our numerical results 
are shown in Fig.[T] Here, we have chosen gt = 106.75 for 
illustration. For most of the parameter range, one can find 



noh' cc — Mg 
(crv) 



-3/2 



(11) 



This feature can be approximately derived from the Boltz- 
mann equation when we omit the equilibrium number density 
Yeq of Eq. ([T) and require p/p^ » 1 at the usual freeze- 
out time of Xf x 20. In Fig. [1] the relic density bound de- 
scribes the standard cosmology case H = Hsc- In this case, 
the relic density bound (crv) > 2 x 10"~''cm^' sec"', indicates 
that the predicted relic density is not bigger than the observed 
value. As shown in Fig. [1] the BBN, unitarity and relic den- 
sity bounds, together, strictly constrain the thermally averaged 
annihilation cross section. 



2 X 10"^* < (crv) < 3 X 10"^-min 



3m TeV^ 
TeV ' rrP- 



(12) 



(P/Pch 



< 1 



(9) where {crv} is in units of cm-' sec 
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B. The quintessence scenario 

In this subsection, we consider another possible scenario 
in which the dark energy is not negligible in the very early 
Universe. It is well-known that in the present Universe, the 
dark energy has begun to play a dominant role. Due to the 
repulsive gravity of dark energy, the Universe is currently un- 
dergoing an accelerated expansion. However, we still do not 
know whether the dark energy is the cosmological constant or 
some dynamical scalar field. If the dark energy is some scalar 
field, for example, the quintessence field, it is possible that 
the dark energy can dominate the early Universe. Of course, 
in the early Universe, the kinetic term dominates the energy 
density of the scalar field. 

In the quintessence scenario, the Friedmann equation can 
be written as isj] 



3^ 



1+^ 



(13) 



where is the quintessence dark energy density. The kinetic 
term domination leads to 



15 10.752 MeV2 



(14) 



The parameter 77 s ('^*'''^'')r-iM v defined by the ratio 
of quintessence-to-photon energy densities at the temperature 
r = 1 MeV. Considering the BBN constraint, we require 



(p<s>/ph 



< 1 . 



(15) 



Then, one can derive the BBN bound for the quintessence sce- 
nario, 



77 < 5.4 



(16) 



If p<t>/p <*c 1 at the usual freeze-out time Xf x 20 (np-rj <s: 
2 X 10 m is in unit of GeV), one can totally omit the dark 
energy contribution. In this case, one can derive the standard 
cosmology bound (crv) ^ 3 x lO^^Scm'' sec"'. However, if 
the quintessence kinetic energy density dominates the early 
Universe, (crv) can be enhanced significantly. Our numerical 
results are shown in Fig.[2l For most of the parameter space, 
there is a relation between the dark matter mass m and annihi- 
lation cross section (crv): 



(crv) log(m X C) 



(17) 



where C ~ y?; x 10^ and m is in units of TeV. This fea- 
ture can also be approximately deduced from the Boltzmann 
equation if one ignores the equilibrium number density Yeq 
and requires m-r] » 2 x 10 One can see that the BBN, 
unitarity and relic density bounds can effectively constrain 
the thermally averaged annihilation cross section {crv}. For 
10 GeV < ;77 < 10 TeV, we obtain 



2 X lO"^*" < (o-v) < 3 X 10""min 



4m 



1 



6 + log(m) ' 

where m is in units of TeV and (crv) is in units of cm-' sec"' 



(18) 




10' 



10 10' 
m (GeV) 



FIG. 2: The allowed region of (crv) and m for different ;/ in 
the quintessence scenario from the observed relic density Cloh' = 
0.1131 ± 0.0034. Here, we also consider the unitarity bound in Eq. 
^ and the BBN bound in Eq. l |16t . The region with dotted bound- 
ary denotes the allowed range of (crv) in Eq. l|20j from the PAMELA 
experiment. 



C. The PAMELA, ATIC and PPB-BETS Anomalies 

Recently, the indirect dark matter detection experiment 
PAMELA lfl2ll reported an excess in the positron fraction from 
10 to 100 GeV, but showed no excess for the antiproton data. 
Dark matter annihilation can account for the PAMELA exper- 
iment if 1 13] 

BF X (crv) ~ 3 X 10"^^^ ^^^'^ ' (^9) 

where BF is the boost factor. Current analysis on the dumpi- 
ness of dark matter structures indicates that the biggest prob- 
able boost factor should be less than 10-20 I14il . Therefore, 
we have 



(o-v) ~ 3 X 10" 



3x10 



-23 / m \2 
I TeV/ 



cm"* sec ' 



(20) 



In the standard cosmology, such a large (crv) will give a 
smaller dark matter abundance than the observed one. So far, 
many authors have provided several very interesting mecha- 
nisms to try to resolve the relic density puzzle, including the 
Sommerfeld enhancement jT^ and the Breit-Wigner enhance- 
ment lfl6ll (for the relevant discussions, see Ref. iIitIi ). In these 
mechanisms, (crv) will increase with the expansion of the Uni- 
verse. If (crv) is a constant, one must resort to the non-thermal 
production of dark matter |fl8[|. In the cuspy halo profile case, 
the annihilation cross section in Eq. ( |20l i will produce abun- 
dant gamma rays from galactic center, which conflicts with 
the HESS results |[l9ll . In fact, current observations do not 
give the direct evidence that any nearby galaxy has a cusped 
dark matter halo profile. In addition, the long-lived interme- 
diate state can relax the gamma ray constraint ||20ll . 
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In the brane-world scenario, the authors of Ref. 112111 ap- 
proximately estimate that the large expansion rate H may 
admit (cry) ~ lO^^GeV"^, which can accommodate the 
PAMELA results and the correct relic density. In this paper, 
we have numerically calculated the Boltzmann equation and 
constrained m, {crv) and M5 (or 77) from the PAMELA results 
in both brane-world and quintessence scenarios. As shown 
in Figs. [T] and |2] (the region with dotted boundary), the re- 
quired (crv) in Eq. (|20] | can produce the observed dark matter 
abundance. It is worthwhile to stress that the unitarity bound 
indicates m < 3 TeV for the PAMELA result. In addition to 
the PAMELA experiment, the ATIC IH] and PPB-BETS IH 
balloon experiments have also seen the excess in the + e' 
energy spectrum between 300 and 800 GeV. In order to ex- 
plain the ATIC and PPB-BETS anomalies, the dark matter 
mass m should be order of TeV. It is clear that the brane- 
world and quintessence scenarios can simultaneously explain 
the dark matter relic density and the PAMELA, ATIC and 
PPB-BETS anomalies. 



III. DISCUSSION AND CONCLUSION 

We have constrained the dark matter thermally aver- 
aged annihilation cross section (crv) in the brane-world and 
quintessence scenarios. In these two scenarios, the Hubble 
expansion rate can be written as 



It should be emphasized that our previous analysis in Sec. II 
can be applied to other models in which Z(r) is proportional 

to r-* or r-. 



In conclusion, in the scenarios of brane-world and 
quintessence, the modified Hubble expansion rate H can en- 
hance the thermal relic density of dark matter. Further- 
more, the constraint on the thermally averaged annihilation 
cross section (crv) ~ 3 x 10"^^ cm"* sec"' can be relaxed 
for the 5-wave annihilation. Considering the relic density, 
BBN and unitarity bounds, we have derived the general con- 
straints on (crv) in Eqs. (fTSl i and (fTsT l for the brane-world 
and quintessence scenarios, respectively. As shown in Figs. 
[1] and |2] the two scenarios discussed in this paper can si- 
multaneously explain the observed dark matter abundance 
fifl/i^ = 0.1131 + 0.0034 and the PAMELA, ATIC, and PPB- 
BETS anomalies. It should be mentioned that the observa- 
tional data from the dark matter may be used to constrain cos- 
mological models. Once the dark matter mass m and anni- 
hilation cross section (crv) are fixed by the dark matter search 
experiments, we can probe the very early stage of the Universe 
from our numerical results in Figs. [T]and|2] 
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